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We present a comprehensive analysis of constraints on the sterile neutrino as a dark matter candi- 
date. The minimal production scenario with a standard thermal history and negligible cosmological 
lepton number is in conflict with conservative radiative decay constraints from the cosmic X-ray 
background in combination with stringent small-scale structure limits from the Lyman-alpha forest. 
We show that entropy release through massive particle decay after production does not alleviate 
these constraints. We further show that radiative decay constraints from local group dwarf galaxies 
are subject to large uncertainties in the dark matter density profile of these systems. Within the 
strongest set of constraints, resonant production of cold sterile neutrino dark matter in non-zero 
lepton number cosmologies remains allowed. 



PACS numbers: 95.35.+d,14.60.Pq,14.60.St,98.65.-r 



I. INTRODUCTION 

The nature of the dark matter remains a fundamen- 
tal problem in cosmology and particle physics. Much 
can be gained from the inferred non-gravitational prop- 
erties of the dark matter such as its decay, annihilation, 
interaction cross-section with baryonic matter, and ki- 
netic properties pj. A class of candidate dark matter 
particles with no standard model interactions, but cou- 
plings to the standard model neutrinos via their mass 
generation mechanism, are the sterile neutrinos. Ster- 
ile neutrinos can be produced in the early universe via 
non-resonant matter-affected oscillations 0, or through 
a resonant mechanism if there exists a non-negligible lep- 
ton asymmetry |3|]. 

Weak interaction singlets such as sterile neutrinos arise 
naturally in most extensions to the standard model of 
particle physics. Grand unified theories commonly con- 
tain singlets, which can act as sterile neutrinos, and can 
play a role in their mass-generation mechanism [i| . Sev- 
eral models contain light singlets as sterile neutrinos, 
including left-right symmetric (mirror) models [f|, su- 
persymmetric axinos as sterile neutrinos |(|, superstring 
models Q , models with large extra dimensions [a El , and 
phenomenological models such as the z/MSM |l0|. Let- 
ting the mass and mixing angle between the sterile and 
respective active neutrino be free parameters, as well as 
the lepton number of the universe be free, the sterile 
neutrino can behave as hot, warm or cold dark matter, 
with masses in the range ~ 0.1 — 100 keV 11]. The 
abundance of sterile neutrinos and its relationship to the 
mixing parameters is affected by the quark-hadron tran- 
sition, however, this relationship is well known for the 
standard prediction of a cross-over transition [T^ | . 

In the standard non-resonant production mechanism, 



sterile neutrinos are produced via a collision-dominated 
oscillation conversion of thermal active neutrinos. De- 
viations from a thermal spectrum in the sterile neutri- 
nos are produced due to the change in the primordial 
plasma's time-temperature relation during production, 
dilution due to the disappearance of degrees of freedom, 
the modification of the neutrino thermal potential from 
the presence of thermal leptons, and the enhanced scat- 
tering rate on quarks above the quark-hadron transi- 
tion |l3( . As such, the resulting dark matter momentum 
distribution is suppressed and distorted from a thermal 
spectrum. 

Sterile neutrinos exhibit a significant primordial veloc- 
ity distribution. This has the effect of damping inhomo- 
geneities on small scales and thus sterile neutrinos behave 
as warm dark matter (WDM). Models with a suppres- 
sion of small scale power have drawn attention due to 
their potential alleviation of several unresolved problems 
in galaxy and small scale structure formation ^4|. Of 
particular interest recently are the possible indications 
of the presence of cores in local group dwarf galaxies, in- 
ferred from the positions of central stellar globular clus- 
ters 0, ^| and radial stellar velocity dispersions • 
The primordial velocity distribution produces a limit to 
the maximum phase-space packing of the dark matter, 
which, if attained, can produce a cored density profile 
for a dark matter halo. Conversely, this places a ro- 
bust limit — the Tremaine-Gunn bound — on the mass and 
phase space of the dark matter particle from observed dy- 
namics in galaxy centers |18| . 

Lighter mass WDM particles more easily escape grav- 
itational potentials, and therefore suppress structure on 
larger scales, which can be constrained by the observed 
clustering on small scales of the Lyman-alpha (Lya) for- 
est. Possibly the most stringent limits on the suppression 
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of power on small scales are placed by inferring the small 
scale linear matter power spectrum from observations of 
the Lya forest [H IH SJ 1H The same flavor-mixing 
mechanism leading to the production of the sterile neu- 
trino in the early universe leads to a radiative decay |23| . 
The decay rate increases as the fifth power of the mass 
eigenstate most closely associated with the sterile neu- 
trino, and increases as the square of the mixing angle, 
producing a lighter mass neutrino and mono-energetic 
X-ray photon at half the dark matter particle mass. X- 
ray observations can either detect or constrain the pres- 
ence of a line flux from surface mass densities of dark 
matter on the sky 0,E3- So far, X-ray observations 
have placed upper limits on the particle mass and mix- 
ing angle relation of the dark matter sterile neutrino with 
observations of the cosmic X-ray background, clusters of 
galaxies, fi eld gala xies, local dwarf galaxies and the Milky 
Way halo HIM El S H3 ■ 

Sterile neutrinos lighter than those that may be the 
dark matter can also play a cosmological role as hot 
dark matter. One or more such light sterile neutrinos 
may be required to produce the flavor transformation 
seen in the Los Alamos Liquid Scintillator Neutrino De- 
tector (LSND) experiment |3lllp|. Such sterile neutri- 
nos would be associated with mass eigenstates of order 
1 eV, and therefore much lighter than a warm or cold 
dark matter sterile neutrino. A light sterile neutrino of 
the type required by LSND would be thermalized in the 
early universe |33l l34l |. and is constrained by limits on 
the presence of hot dark matter from measures of large 
scale structure [35l |3(|. Such limits can be avoided if 
the LSND-type sterile neutrino was not thermalized due 
to the existence of a small lepton number, though they 
nonetheless may be produced resonantly |37j . 

There are two other interesting physical effects when 
sterile neutrinos have parameters such that they are 
created as the dark matter in the non-resonant pro- 
duction mechanism. First, asymmetric sterile neutrino 
emission from a supernova core can assist in producing 
the observed large pulsar velocities above 1000 km s _1 
jssl Isfll Ecj . The parameter space overlaps that of the 
non-resonant production mechanism (Fig. Second, 
the slow radiative decay of the sterile neutrino dark mat- 
ter in the standard production mechanism can augment 
the ionization fraction of the primordial gas at high- 
redshift (high-z) |4l|. This can lead to an enhancement 
of molecular hydrogen formation and star formation, but 
also results in a strong increase in the temperature of 
the primordial gas |42j . This effect may have dire conse- 
quences on the formation of the first stars, which remains 
an open question |43[ . 

In this paper, we review all of the constraints on the 
parameter space of the sterile neutrino as a dark matter 
candidate, in conjunction with the parameters needed for 
oscillation-based resonant and non-resonant production 
mechanisms for the sterile neutrino as the dark matter. 
In jJETl we review the best current constraints on the ster- 
ile neutrino as a dark matter candidate. We review X- 



ray observation constraints in £111 Al The most promising 
upper mass constraints come from X-ray observations of 
local group dwarfs, but we show that they are subject 
to significant uncertainties in the dark matter profile of 
the dwarf. In £ 111 Bl we discuss constraints from observa- 
tions of the Lya forest. This provides the most stringent 
lower mass constraints, which requires the free streaming 
length to be < 90 kpc/h, and forces the primordial veloc- 
ity dispersion to be too small to sustain cored dark mat- 
ter halos. With the combined constraints, we show that 
the standard zero-lepton number non-resonant produc- 
tion model is excluded if the most stringent constraints 
from the Lya forest are combined with the most conser- 
vative decay limits of the X-ray background, and cannot 
be evaded in a model that dilutes and cools the dark 
matter sterile neutrino by massive particle decay |44| . 
However, we show that the combined constraints do not 
exclude resonant production mechanisms. 



II. CONSTRAINTS 



The non-resonant "zero" lepton- number production 
calculation presented in Ref. [Lj is the simplest case 
model for the production of sterile neutrinos as dark 
matter candidates. In this model, there are no extra 
couplings postulated for the sterile neutrinos, and the 
cosmological lepton number is comparable to the baryon 
number, and thus negligible. The thermal history (up 
to temperatures of T ~ 500 MeV) is given by lattice 
QCD calculations through the quark-hadron transition, 
and contains no sterile neutrinos in the initial conditions 
of the plasma, due to the fact that the thermal potential 
suppresses sterile production at high temperatures. 

In Fig. ^ we show contours of constant comoving den- 
sity comparable to the dark matter density for the non- 
resonant zero lepton number (L = 0) well as 
enhanced resonant production with initial cosmological 
lepton number cases (L = 0.003, L = 0.01, L = 0.1) 
from Ref. We have labeled the standard prediction 
of L ~ 10~ 10 as nil since it is negligible for the non- 
resonant production mechanism. We define the lepton 
number as 



L = 



(1) 



where n v (rip) is the number density of the neutrino (an- 
tineutrino) flavor with which the sterile is mixed, and n 7 
is the cosmological photon number density. The cosmo- 
logical lepton number is limited by the inferred primor- 
dial helium abundance and the large to maximal mix- 
ing angle solutions to the solar and atmospheric neutrino 
problems: \L e \ < 0.05, \L„ + L T \ < 0.4 [H E3 Eg. 
The constraints discussed below are framed around the 
parameter space required for the sterile neutrino dark 
matter production, and many are shown in Fig. ^ 
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FIG. 1: Full parameter space constraints for the sterile neu- 
trino production models, assuming sterile neutrinos consti- 
tute the dark matter. Contours labeled with lepton number 
L — 0, L = 0.003, L = 0.01, L = 0.1 are production predic- 
tions for constant comoving density of Q a = 0.24 for L = 0, 
and Q s — 0.3 for non-zero L |l2||. Constraints from X-ray ob- 
servations include the diffuse X-ray background (green) [27) . 
from XMM-Newton observations of the Coma and Virgo clus- 
ters (light blue) |28| . The diagonal wide-hatched region is the 
claimed potential constraint from XMM-Newton observations 
of the LMC 0. The region at m s < 0.4 keV is ruled out by 
a conservative application of the Tremaine-Gunn bound fl4| . 
The regions labeled Lya are those from the amplitude and 
slope of matter power spectrum inferred from the SDSS Lya 
forest [ Lya (1)] |21|. using high-resolution Lya data [Lya 
(2)] Hi [Si, and that from the high-z SDSS Lya of SMMT 
[Lya (3)] |22|. The grey region to the right of the L = case 
is where sterile neutrino dark matter is overproduced. Also 
shown is the horizontal band of the mass scale consistent with 
producing a 100 - 300 pc core in the Fornax dwarf galaxy [45) . 
The parameters consistent with pulsar kick generation are in 
horizontal hatching [38l l39l |40| . 



A. X-ray measurements 

In this section, we review the sterile neutrino dark mat- 
ter constraints that come from measurements of the X- 
ray background, X-ray measurements from the Virgo and 
Coma clusters, as well as measurements of X-ray fluxes 
from the Draco local group dwarf. 

It is straightforward to translate X-ray astronomy mass 
and mixing angle parameter space constraints to sterile 
neutrino mass constraints in the simplest model with the 
inversion of the production relation in Ref. |13| : 



sin 2 20 = 7.31 x lfW-^ 
VkeV/ 



-1.63 



0.26 



0.813 



(2) 



where 9 is the mixing angle between the active and sterile 
flavor states, and fl s is the fraction of the cosmological 
critical density in sterile neutrinos. 

If sterile neutrinos constitute the density associated 
with dark matter, then their radiative decay would lead 
to a contribution to the X-ray background (XRB here- 
after) [2^|. In a recent analysis of the observed X-ray 
background from HEAO-1 and XMM-Newton, Boyarsky 
et al. 0| place the following limit on the particle mass 
and mixing angle, 



sin 2 29 < 1.15 x 10 



_ 4 ( m s V 
VkeV/ 



0.26 



(3) 



with the corresponding exclusion region shown in Fig. ^ 
If we combine this result with the production Eq. @ , we 
find the corresponding upper mass limit to be 



m s < 8.9 keV. 



(4) 



More stringent limits can be placed by X-ray obser- 
vations of the large dark matter surface mass density 
in clusters of galaxies [2(j. A recent analysis of XMM- 
Newton observations of the Virg o and Coma clusters was 
presented in Boyarsky et al. [28| . More specifically, it was 
shown that near the center of the Virgo cluster (at radial 
distances r < 11 arcmin), the X-ray flux places a rough 
power-law constraint on the sin 2 29 — m s plane, as 



sin 2 20 < l0- 2 ( T " 



' k. v) 



-6.64 



(5) 



If this result is combined with the production mass- 
mixing angle relation [Eq. @], it results in a particle 
mass constraint of |6J| 



< 10.6 keV. 



(G) 



The combined Virgo and Coma analysis of Boyarsky 
et al. [2^ presents a more stringent limit. In this case, 
an approximate power-law fit to their exclusion region 
places a limit of 



sin 2 29 < 8 x 10~ 5 ( 

VkeV/ 



-5.43 



(7) 



which is considerably stronger than the XRB limit, 
Eq. ©. Using the production relation Eq. J2J), this limit 
yields a sterile neutrino mass limit of 



< 6.3 keV, 



(8) 



an improvement on the XRB limit, Eq. (f4"|). 

There are some notable issues with the Boyarsky et 
al. j2^| analysis of the Virgo and Coma cluster data. The 
analysis uses a fixed phcnomcnological model for the X- 
ray emission of the cluster with specific lines added to fit 
atomic lines in the spectrum. On top of the phenomeno- 
logical model, a Gaussian line representing the potential 
sterile neutrino flux is inserted with the width of the en- 
ergy resolution of the instrument. This method does not 
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allow modeling of the energy-dependent resolution of the 
response of the detector, and can lead to a non-detection 
of a line feature that exists at the position of an atomic 
line. Emission lines could be more properly modeled 
for the gas in clusters by using a Mewe-Kaastra-Liedahl 
(MEKAL) model of the atomic and bremsstrahlung emis- 
sion of the gas [4!j. Nevertheless, barring the chance 
coincidence of the sterile neutrino emission feature lying 
on an instrumental feature or an atomic line, the limits, 
Eq. H7I8|) from the Coma plus Virgo analysis of Ref. 
should be robust. 

We now discuss the prospects of detecting sterile neu- 
trino dark matter in local dwarf galaxies. It was re- 
cently proposed by Boyarsky et al. [2j| that X-ray ob- 
servations by XMM-Newton of nearby local group dwarf 
galaxies may present the best opportunity for constrain- 
ing or detecting the sterile neutrino decay flux of X-ray 
photons. The constraint region from that paper, us- 
ing XMM-Newton observations of the Large Magellanic 
Cloud (LMC), is shown as the broad diagonal hatched 
region in Fig. ^ 

Any analysis of the expected X-ray flux from sterile 
neutrino decays, such as in Ref. [2^, is prone to the un- 
certainties in the dynamical estimate of the dark matter 
distribution in the dwarf galaxy [Hcj j. For example, in 
the analysis of Ref. the dark matter profile of the 
Draco dwarf galaxy is consistent with both cored and 
cusped dark matter distributions such as the NFW |§2 
and Burkert |53|| profiles. These profiles may arise in 
the case of massive "cold" sterile neutrinos (NFW) or 
for lighter "warm" sterile neutrinos. In order to demon- 
strate the uncertainties in the X-ray flux due to the dark 
matter distribution in Draco, we show in Fig.[21the value 
of the quantity J[Af2(0)] which is defined as the line of 
sight integral of the matter distribution over a solid angle 
Afi(0) centered on the dwarf galaxy and expressed as 

j [An(0)] = 











/ * 




/ I[f{x)} dx 


d9'. (9) 


10 


Jo 


^min(e') 





Here, /[f (a:)] is a function which depends on the assumed 
dark matter profile, and takes the form of 
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NFW 



'BUR 



[r»] 



[r»] = 



1 



f{x) [1 + f(x)] 
1 

[l + r(x)] [l + f 2 (x)} 



(10) 
(11) 



for the NFW and Burkert profiles respectively. In Eq. 10, 
p s is the characteristic density of the assumed profile, 
the limits of the integration along the line of sight are 
£max,min# = D cos 6* ± r 2 + (D sine) 2 , the quantity f is 
defined as f(x) = r(x)/r s , where r s is the scale radius 
of the assumed profile, and it relates to the line of sight 
element through f(x) = \J x 2 + D 2 — 2xD cos e. The dis- 
tance D to Draco is D — [75.8 ± 7 ± 5.4] kpc estimated 
using RR Lyrae variable stars [54| , and the tidal radius 
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FIG. 2: The value of the quantity J[AQ.(6)] (see text) for two 
representative dark matter distributions in the Draco dwarf 
galaxy. The thickness of both curves corresponds to the range 
of values in each profile due to the distance uncertainties to 
Draco. The top curve corresponds to an NFW profile with 
log(p s /M Q kpc" 3 ) = 7.0 and log(r s /kpc) = 0.85, while the 
lower curve depicts a Burkert profile with log(p s /M0kpc -3 ) = 
9.0 and log(r 3 /kpc) = -0.75. 



r t is taken to be r t 7 kpc. However, we point out that 
as long as r t ^> r s the dependence of Afl(e) on r t is weak. 

In Fig. |21 we show the value of the quantity J[Af2(6*)] 
for two representative cases of a cored and cusped dark 
matter distributions. The top band depicts an NFW pro- 
file with log( j o s /Af Q kpc -3 ) = 7.0 and log(r s /kpc) = 0.85, 
while the bottom band corresponds to a Burkert profile 
with log(p s /M Q kpc -3 ) = 9.0 and log(r s /kpc) = -0.75. 
Both of these profiles are consistent with the observed 
line-of-sight velocity dispersion measurements as shown 
in |5l) . Based only on profile and distance uncertain- 
ties, the value of the quantity Af2(0) varies between 
[2.5 x 10 17 - 1.3 x 10 18 ]GeVcm- 2 , a factor of 5.2 at 8 
arcmin (which corresponds roughly to the XMM field of 
view) . 

It has been claimed by Boyarsky et al. 29] that per- 
haps the X-ray flux from the LMC could provide the 
strongest constraint on the parameters of a sterile neu- 
trino dark matter candidate. Modeling the distribution 
of dark matter in the inner regions of the LMC is even 
more uncertain due to the presence of a stellar disk and 
a bar (as for the Milky Way as well). As it was shown 
in numerous studies, e.g., Ref. |55j . the LMC is baryon 
dominated in the central region, with a mass-to-light ra- 
tio of ~ 3 within the inner ~ 9 kpc (for comparison, 
dwarf spheroidals have mass-to-light ratios of ~ 100). 
The importance of baryon domination on the distribution 
of dark matter in the LMC was shown in the analysis of 
Ref. j5|J . The derived mass of the LMC is uncertain to 
within 20%, depending on whether the disk is modeled as 
"maximal," or "minimal." If the baryons are dominat- 
ing, then the rotation curve is much less sensitive to the 
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distribution of dark matter, making any estimate of the 
dark matter mass of the the LMC unreliable. In light of 
these uncertainties, we conclude that the LMC is an un- 
reliable Milky Way satellite for robust X-ray constraints 
from detection of sterile neutrino decays. 

Blank sky observations by X-ray telescopes may also 
provide a detectable dark matter decay flux due to the 
dark matter halo of the Milky Way itself [H H3- The 
importance of uncertainties in the dark matter profile 
of the Milky Way halo can also be significant. Models 
of the measures of the dynamics of the Milky Way can 
fit a range of halo masses of (0.7 - 2) x 10 12 M Q [EJ 
l5^ |. The uncertainty in these dynamical estimates can 
lead to a factor of 3 difference in the expected X-ray flux 
in directions perpendicular to the galactic plane, though 
these uncertainties are not reflected in the analyses of 
Refs. HHU. 

B. Ly man-alpha Forest 

The standard paradigm of cosmological structure for- 
mation is the gravitational growth and eventual collapse 
from small to large scales of initially adiabatic, Gaussian 
density fluctuations. WDM particle candidates alter the 
initial conditions of the perturbation spectrum by damp- 
ing small scale fluctuations below the free streaming scale 
of the WDM. The sterile neutrino particle mass is con- 
strained from below by the observations of small scale 
cosmological structure. 

The most stringent lower bounds on the sterile neu- 
trino mass arise from observations of the clustering of 
gas along the line of sight to distant quasars. The den- 
sity fluctuations of the gas follow that of the dark matter 
to the scale at which the gas becomes pressure supported. 
The density fluctuations are linear to mildly nonlinear, 
and can probe extremely small scale dark matter fluctu- 
ations. 

Using a combination of cosmic microwave background 
observations, the shape of the 3D power spectrum of 
galaxies from the Sloan Digital Sky Survey (SDSS), the 
inferred linear matter power spectrum amplitude and 
slope from the SDSS Lya forest, the lower limit on the 
sterile neutrino particle mass is m s > 1.7 keV (95% CL). 
In this limit, a linear bias relation between the Lya 
flux power spectrum and matter power spectrum was 
assumed |2l|. Using the inferred matter power spec- 
trum from high- resolution spectra of the Lya forest, the 
limit was improved to m s > 3.0 keV (95% CL); how- 
ever, significant systematic uncertainties exist in mod- 
eling the high-resolution data |2l|. Seljak et al. [2^ 
(hereafter SMMT) find a much more stringent constraint 
when directly using high-z flux power spectra from the 
SDSS measured by McDonald et al. [5!| and other higher- 
resolution flux power spectra: 

m s > 14 keV (95% CL). (12) 

All of these Lya constraints are shown in Fig. ^ 



Pioneering work on the sterile neutrino dark matter 
transfer function, using the approximation of a sup- 
pressed thermal sterile neutrino distribution, was done 
by Refs. HJ|6l|. It has been shown that the sterile neu- 
trino momentum distribution is not simply a suppressed 
thermal distribution, but is significantly nonthermal due 
to the effects on production of the dark matter by the 
changing particle population in the early universe, lep- 
ton population affecting the neutrino thermal potential, 
the quark-hadron transition, and the dilution of the dark 
matter due to particle annihilation 

GSm. A11 of the 

above Lya constraints use either the appropriate non- 
thermal transfer function for sterile neutrinos modified 
by their production at high temperatures in the early uni- 
verse at or near the quark-hadron transition 0] ■ SMMT 
include an approximation to these effects via an augmen- 
tation in the sterile neutrino particle mass of 10%, since 
the above effects cool the momentum distribution. The 
original assumption of a simple suppressed thermal dis- 
tribution for the sterile neutrino produces a suppression 
scale that is altered by a factor of order 10% in the par- 
ticle mass of the sterile neutrino, though the correction 
increases with sterile neutrino particle mass due to their 
production at higher temperatures where all of the above 
effects are more pronounced. All of the physical effects 
producing a nonthermal distribution where included in 
the transfer function fit given in Ref. |l3j . 

The strength of the SMMT results arises from the sen- 
sitivity of the high- z Lya flux power spectra to changes 
in the ID linear matter power spectrum, which is itself 
more sensitive to the suppression scale of WDM than the 
3D power spectrum. At high-z, the recovery of the am- 
plitude of the power spectrum via nonlinear clustering is 
reduced, enhancing the effects of WDM suppression. The 
temperature-density relation of the gas is constrained si- 
multaneously by the observations in SMMT, though the 
strong change in the thermal state of the gas due to radia- 
tive sterile neutrino decay may be significant 42]. An- 
other essential feature of the analysis in SMMT is the 
use of smaller volume hydrodynamical simulations that 
can resolve the very small free-streaming scale of 14 keV 
neutrinos. The free streaming scale of a sterile neutrino 

WDM is [m 

where m s is the mass state closely associated with 
the sterile neutrino flavor, and (p/T) is the mean mo- 
mentum over temperature of the sterile neutrino dis- 
tribution, and h is the Hubble parameter in units of 
100 km s _1 Mpc -1 . A thermal WDM particle is de- 
fined such that (p/T)/3.15 s» 1. Due to the thermal 
history of the universe during production, the standard 
non-resonant production mechanism produces a "cool" 
sterile neutrino distribution, (p/T)/3.15 s» 0.9 while 
the resonant production mechanism enhances low-p pro- 
duction, and (p/T) /3. 15 « 0.6, depending sensitively on 
the mass of the neutrino and initial lepton number [111 ]. 
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The mass scale within the free streaming length is 

'°Mw)(^)W 

where Q m is the fraction of the cosmological critical den- 
sity in matter. 

To resolve the required Aps ~ 54 kpc//i of a 14 keV 
neutrino, SMMT use a 20 Mpc/i _1 box with 256 3 par- 
ticles in dark matter and 512 3 cells for gas, providing 
a grid spacing of 39kpcft.~ 1 . This allows only for a 
resolution of a fraction of the suppression due to free 
streaming. Higher resolution in principle should only en- 
hance the effects of the WDM suppression. To test con- 
vergence, SMMT use a single smaller volume simulation 
(10 Mpc/i" 1 ) with the same particle and grid spacing 
and find a ~20% change in the magnitude of the effect, 
though no higher resolution simulations were performed 
to test if the change is subsequently smaller as would 
be expected in numerical convergence. In addition, two 
other physical effects are degenerate, to different extents, 
with the effects of WDM, namely, pressure support due to 
the Jeans scale of the gas, and temperature broadening. 
These can mimic or hide the effects of the free streaming 
scale of the WDM, therefore the thermal properties of 
the gas are crucial to properly model in such limits. 

A recent analysis by Viel et al. 62] (herafter VLHMR) 
used the same data of the SDSS Lya flux power spectrum 
of McDonald et al. j5(| , but excluding higher resolution 
data used by SMMT. VLHMR utilizes a different method 
of mapping the response of the flux power spectrum to 
changes in cosmological and astrophysical parameters, 
viz., VLHMR uses a parameterized Taylor series expan- 
sion of the flux power response to changes in physical 
paramters, with the Taylor parameters fit by numerical 
simulations. This method was shown to give similar re- 
sults for the standard ACDM cosmology ]63| . VLHMR 
finds a weaker sterile neutino particle mass lower limit, 
10 keV (95% CL), than that of SMMT, 12 keV (95% CL), 
using the Lya data from SDSS alone. The discrepancy 
with SMMT is increased when taking into account that 
the limit from VLHMR is too strong by >10% because 
they do not use the correction for the "cooler" nonther- 
mal spectrum of the sterile neutrinos due to the effects 
in the early universe described above. Overall, between 
the two analyses there exists a discrepancy of > 30% 
in the limits on the sterile neutrino mass. This could 
arise due to different CMB and galaxy data sets used: 
SMMT employs WMAP first year data and the SDSS 
galaxy power spectrum; VLHMR employs WMAP third 
year data, higher resolution CMB measurements and the 
2dF galaxy power spectrum. Therefore, it is not certain 
whether the different CMB and galaxy data sets are the 
source of the discrepancy, or whether it is in the hydrody- 
namical simulations and method of mapping the response 
of the flux power spectrum to changes in physical param- 
eters. However, both analyses find a stringent limit due 
to the precision of the measurement of the McDonald et 



al. |59j flux power spectrum at high redshift. 

To reflect back on one of the principle motivations for 
WDM, it is important to note that it was show n by Stri- 
gari et al. ^5| that if the constraints from Ref. (2J , and 
especially SMMT are valid, then dynamical constraints 
from the Fornax dwarf galaxies limit the size of a core 
to be < 85 pc and -C 10 pc for the two Lya limits, re- 
spectively. In Fig. n we show the particle mass required 
to produce a [100-300] pc core b y in verting the dynam- 
ical constraints of Strigari et al. 45]. These constraints 
by Strigari et al. are based on the inferred dark matter 
density profile from the radial velocity dispersion profile. 
The positions of the globular clusters in Fornax may indi- 
cate a core of -240 pc EH EH- Furthermore, the SMMT 
particle mass limit also limits the scale of the suppression 
of the halo mass scale to be well below the typical masses 
of dwarf galaxies, Mps < 10 7 M Q h _1 . 

In the resonant production model, the exact level of the 
lower bound from the Lya forest would be modified for 
each lepton number case due to variation of (p/T)/3.15 
and therefore Afs for each case. However, this is at the 
level of —30% and is not monotonic across the lepton 
number region in m s — sin 2 29 space. Therefore, we leave 
Lya forest as a horizontal line in Fig. ^ to provide a 
rough guide to the limit. As shown in Fig.QJ even with 
the inclusion of all constraints, the resonant production 
model remains unconstrained at high-mass scales: 14 < 
m s < 100 keV. 

To summarize, the results of SMMT are extremely sig- 
nificant, ruling out much of the parameter space that 
motivates WDM in general, and when combined with 
conservative X-ray bounds, as shown above, they rule 
out the standard sterile neutrino production mechanism. 
Therefore, there is strong motivation to verify the robust- 
ness of the SDSS Lya measurements employed by SMMT 
as well as the modeling by hydrodynamic simulations. 



III. CONSTRAINTS IN A DILUTION 
SCENARIO 

It has been proposed that the production of sterile neu- 
trino dark matter could be followed by the decay of a 
massive particle, whose decay products reheat the cou- 
pled species in the plasma, dilute the sterile neutrino dark 
matter and cool it relative to the coupled species [44| . 
Though this involves a conspiracy between the lifetime 
of the massive species and parameters coupling the sterile 
neutrino to the active sector, it is an interesting possibil- 
ity that may alleviate structure formation constraints on 
the sterile neutrino. However, it is important to note 
that it does not allow a window for warm dark matter, 
as this mechanism cools the WDM particle until it may 
be consistent with structure formation limits. 

In this scenario, a massive particle (m — 100 GeV) 
decays so that the entropy release changes the relative 
abundance of the dark matter by a factor S, i.e. fi s — > 
fl s /S. In the non- resonant oscillation production model, 



7 



100 



Co 10 



Lya 



X— ray Background 



10 

m s [keV] 



15 



20 



TO c„tropy > 3 ke y^ ^g) 

This is in conflict with the XRB in the entropy release 
model, Eq. (|17fl . The scaling of these relations is illus- 
trated in Fig. [3J More stringent limits from X-ray clus- 
ters or local group galaxies would be in stronger conflict 
with Eq. (|19() . Blank-sky observations including contri- 
butions from the Milky Way halo also constrain high S 
models |30|. As such, we conclude that massive particle 
decay does not open a new window for the non-resonant 
oscillation production scenario. 



IV. CONCLUSIONS 



FIG. 3: Shown here are the constraints on the massive particle 
decay dilution model. The diagonally-hatched (blue) region 
is the lower-mass Lya limit of Ref. 21], while the vertically 
(red) hatched region the Lya limit of SMMT. In combination 
with the conservative XRB limit (green) [27|. even extreme 
dilution models of S — 100 are in conflict with combined 
constraints. The standard case of no dilution corresponds to 
S=l. 

requiring a subsequent increase in the mixing angle of the 
production relation to 

/ m N -1.63 / O \ 813 

(15) 

from Eq. Using this and the conservative XRB limit, 
Eq. J3J), the limit on the entropy release factor is 

^ / /-v \ — 2.23 

S < 8 - 6 * 10 W (55) ■ < K > 

which is valid for the mass range 1 < m s < 100 keV. 
With extended models, S can be as large as 100, at which 
point the dilution is occurring for sterile neutrinos that 
were at or nearly at thermal equilibrium with the plasma 
prior to the massive particle decay. Using Eq. (|16fl and 
the case where S — 100, the limit on sterile neutrino dark 
matter from the XRB is 

TO o„tropy < 2 g ke y ( 17 ) 

Limits from cosmological structure such as the Lya 
forest are also modified by the cooling of the sterile neu- 
trino free streaming length due to entropy release, such 
that the new lower mass limit is 

TO o„tropy = TO sta„dard £-1/3 ^gj 

For the most stringent Lya forest limit of SMMT, 
Eq. i|12|) . the entropy release model lower limit at S = 100 



We have provided a comprehensive analysis of con- 
straints on the parameter space of interest for the non- 
resonant and resonant production mechanisms of sterile 
neutrino dark matter. Observations in the X-ray of clus- 
ters of galaxies and the XRB place a limit to the radiative 
decay rate of a sterile neutrino candidate and provide 
constraints in the upper mass scale of the sterile neu- 
trino dark matter. We have shown that limits from local 
group dwarf galaxies are subject to large uncertainties 
in the dark matter profile of these objects. The lack of 
the effects of the suppression of small scale power in the 
Lya forest stringently limits the low mass scale region of 
parameter space. 

The SDSS Lya constraints from the analysis of SMMT, 
when combined with X-ray constraints, are in conflict 
with the sterile neutrino being the dark matter in the 
standard non-resonant zero lepton number production 
model. The SMMT limits also exclude much of the pa- 
rameter space of interest for general WDM models. Non- 
zero lepton number cosmologies remain allowed for res- 
onant production of "cold" sterile neutrino dark matter. 
We find, however, that dilution scenarios do not open a 
window for sterile neutrino dark matter. If the X-ray and 
Lya constraints remain robust, then only non-zero lep- 
ton number cosmologies remain viable for the oscillation- 
production models of sterile neutrino dark matter. 
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